Formation of deeply bound molecules via chainwise adiabatic passage 
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We suggest and analyze a novel technique for efficient and robust creation of dense ultracold 
molecular ensembles in their ground rovibrational state. In our approach a molecule is brought to the 
ground state through a series of intermediate vibrational states via a multistate chainwise Stimulated 
Raman Adiabatic Passage (c-STIRAP) technique. We study the influence of the intermediate states 
decay on the transfer process and suggest an approach that minimizes the population of these states, 
resulting in a maximal transfer efficiency. As an example, we analyze the formation of 87 Rb2 starting 
from an initial Feshbach molecular state and taking into account major decay mechanisms due to 
inelastic atom-molecule and molecule-molecule collisions. Numerical analysis suggests a transfer 
efficiency > 90%, even in the presence of strong collisional relaxation as are present in a high 
density atomic gas. 
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Ultracold molecular gases open possibilities for stu- 
dyng new exciting physical phenomena and their appli- 
cations. For example, ultracold molecules can find use in 
testing fundamental symmetries 0, 0] , in precision spec- 
troscopy [1, IU and ultracold chemistry Q . Dipolar ultra- 
cold quantum gases promise to show new phenomena due 
to strong anisotropic dipole-dipole interactions. Dipolar 
molecules in optical lattices can be employed as quantum 
simulators of condensed matter systems @. Ultracold 
polar molecules also represent an attractive platform for 
quantum computation [7[. 

Dense samples of molecules in their ground rovibra- 
tional state v — 0, J = are required for many of these 
applications. In this state, they have a large permanent 
electric dipole moment and are stable with respect to col- 
lisions and spontaneous emission. Currently translation- 
ally ultracold (100 nK - 1 mK) molecules are produced 
by magneto- [H and photoassociation 0] techniques. In 
both of these techniques the molecules are translationally 
cold, but vibrationally hot, since they are formed in high 
vibrational states near the dissociation limit of the elec- 
tronic ground state. Therefore, once created, molecules 
have to be rapidly transfered to the ground rovibrational 
state. 

One of the most efficient ways to transfer population 
between two states is based on the Stimulated Raman 
Adiabatic Passage (STIRAP) technique [H HO, Q3, • 
STIRAP provides a lossless robust transfer between an 
initial and a final state of a three-level system using a Ra- 
man transition with two counterintuitively ordered laser 
pulses. The main difficulty with a two-pulse STIRAP 
in molecules is to find an intermediate vibrational state 
of the excited electronic potential with a good Franck- 
Condon overlap with both a highly delocalized initial 
high vibrational state and a tightlylocalized v = state 
[14] . It was therefore proposed in 15j to transfer popula- 
tion in several steps down the ladder of vibrational states 
using a sequence of stimulated optical Raman transitions. 



In this case the initial and final vibrational levels of each 
step do not differ significantly, and it is easier to find 
a suitable intermediate vibrational level in the excited 
electronic state. In this step-wise approach population 
is transferred through a number of vibrational levels in 
the ground electronic state. In a dense gas, molecules in 
such states are subject to inelastic collisions with back- 
ground atoms or other molecules. The released kinetic 
energy greatly exceeds the trap depth resulting in loss of 
both molecules and atoms from the trap. This process 
is expected to limit the efficiency of creation of dense 
ultracold molecular samples. In this work we present 
a technique allowing an efficient transfer of a molecule 
from a high-lying to the ground vibrational state which 
minimizes population loss due to inelastic collisions in in- 
termediate levels. Our technique is based on generalized 
chainwise STIRAP, which in principle allows for lossless 
transfer to the ground vibrational state. We note that 
serial STIRAP as in [l3l . [Tol and pump-dump technique 
with a train of short pulses [16[ should also allow lossless 
transfer if pulses are shorter than the collisional relax- 
ation time. 

The idea of this work can be described using a sim- 
ple five-level model molecular system with states chain- 
wise coupled by optical fields as illustrated in Fig.[T] 
The states \gi), \gz) and \g%) are vibrational levels of 
the ground electronic molecular state, while \e\) and \e2) 
are vibrational states of an excited electronic molecular 
state. Molecules are formed in a high vibrational state 
| <7i) , which in the following ia assumed to be a molecu- 
lar Feshbach state. The state \ga) is the deepest bound 
vibrational state v = 0, and I02) is an intermediate vibra- 
tional state. The goal is to efficiently transfer population 
from the state |<7i) to state \gs). At least two vibrational 
levels I ei) and |e2) in an excited electronic state are re- 
quired, one having a good Franck-Condon overlap with 
I gs), and the other with the initial Feshbach molecular 
state \gi). 
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FIG. 1: Schematic showing the multistate chainwise STIRAP 
transfer of population from the Feshbach |<?i) to the ground 
1 173) vibrational state. 



In the states |ei ) and |ea) s molecules decay due to spon- 
taneous emission and collisions, and in the states \g±) (for 
bosonic molecules) and |<? 2 ) they experience fast inelas- 
tic collisions with background atoms leading to loss of 
molecules from a trap. It means that populating the 
states \e\), \e 2 ) and | (72) has to be avoided when a back- 
ground atomic gas is present, or the transfer process has 
to be faster than the collisional relaxation time. 

This can be achieved via chainwise STIRAP. The wave 
function of the system is j*) = J2i Q ex P ( — N)> 
where i = g 1 ,e 1 ,g 2 ,e 2 ,g 3 ; <P 9l = 0, (f> ei = v\t, <t> 92 = 

(V2-Vl)t, (f) e2 = (V3+V2-Vl)t, (f>g 3 = {vi ~ V Z + V 2 ~ V X )t \ 

Vi is the frequency of the ith optical field. The evolution 
is then governed by the Schrodinger equation 



91*) 

ih-L± = H(t) 1*) 



(1) 



where the time-dependent Hamiltonian is given by 



H = 
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Here £l;(i) = Pi£i(t)/2h, i = 1,2,3,4 are the Rabi fre- 
quencies of optical fields; is the amplitude of ith optical 
field, [ii is the dipole matrix element along the respective 
transition, Ai = uj\ — V\ and A2 = 0^4 — ^4 are one-photon 
detunings of the fields, and the u>i are the molecular fre- 
quencies along transition i. We assumed in Eq. ([2]) that 
pairs of fields coupling two neighboring ground state vi- 
brational levels are in a two-photon (Raman) resonance. 

The Hamiltonian Eq. |(3J) has a dark state, a specific su- 
perposition of states uncoupled from applied laser fields, 
given by the expression 



^2^4 \gi) - ^4^1 152) + ^1^3 133) 



(3) 



In c-STIRAP (as in classical STIRAP) the optical fields 
are applied in a counterintuitive way, i.e. at t — —00 only 



a combination of the Hi, D3, f2 2 fields, and at t = +00 
only of 1^3, 2 , and f2i is present. As a result the dark 
state is initially associated with the \gi) and finally with 
the |<?3) state. Adiabatically changing the Rabi frequen- 
cies of the optical fields so that the system stays in the 
dark state during evolution, one can transfer the system 
from the initial high-lying \gx) to the ground vibrational 
^3) state with unit efficiency, defined as the population 
of the I33) state at t — +00. The dark state does not have 
contributions from the |ei) and |e 2 ) excited states, which 
means that they are not populated during the transfer 
process. As a result, the decay from these states does 
not affect the transfer efficiency. To minimize the pop- 
ulation in the intermediate ground vibrationl state we 
apply the fields such that SI2, ^3 ^S> fii, O4 and JTq, ^3 
temporally overlap both the f^i and pulses 17]. In 
this case Eq.Q indicates that the population in the \g 2 ) 
state can, in principle, be completely suppressed at all 
times. This is the main idea of this work. 

We now analyze this system in detail. To simplify the 
analysis, we assume Q 2 = ^3 = ^0; ^0 is independent 
of time (in practice the corresponding pulses just have 
to be much longer than f2i(i), fl^t) and overlap both 
of them), and fin S> | ST2 4 1 . We also set one-photon 
detunings to zero Ai = A 2 = 0, and define the effective 
Rabi frequency Q(t) = yfi^+Tq and a rotation angle by 
tan#(t) = rii(t)/Q,4,(t). The eigenvalues of the system (2) 
are £n = 0, corresponding to the dark state, and — 
±f2/i/2, and £3,4 = ±\/2Qq to bright states. Adiabatic 
eigenstates |$) = {|$")}, n = 0, ...4 and the bare states 
are transformed as \&j = ^ ■ Wij$j via a rotation matrix 
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and 



where £ = f2/fin and terms of the order of 0(£ 2 ) 
higher are neglected. 

The adiabaticity condition in this case requires 9 <C 
il, f2n. If the condition 9 <C fi is satisfied, the dark 
state will not couple to the |<j> 1,2 ) states, corresponding 
to the closest in energy £i i2 eigenvalues. Coupling to the 
|$ 3 ' 4 ) states will be suppressed even more strongly, since 
£1 <C r2o- This gives a standard STIRAP adiabaticity re- 
quirement QTtr 3> 1, where T tr is the c-STIRAP transfer- 
time. 

To study the effect of the decay from \g 2 ) and \gi) on 
the dark state evolution, we turn to a density matrix 
description [T^], and use the adiabatic basis states. The 
density matrix equation then takes a form 



do a 



W T W,p a 



na a 

- C p , 



(•5) 



where the density matrix p and the Liouville operator 
C in this basis are given by p a = W T pW and C a p a = 
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W T £pW, and the Hamiltonian H a is diagonal; W is the 
rotation matrix. The Liouvillc operator C consists of the 
usual decays, where only population decays (oc Ty ) into 
other vibrational states or the continuum are considered 
(see FigfT]). Since at t = — oo all population is assumed 
to be in state |<7i), initial conditions for Eq. ([5|) read as 
Poo = 1) Pnm = for nm ^ 00, where pgo denotes the 
dark state population. 

The decay of the dark state due to the population loss 
from the l^) and \g 2 ) states is then described by the 
equation (keeping only terms up to the f2 2 /ri§ order) 

% w -(r 2 + r lC os 2 0) (JL S m20) -r lC os 2 0. (6) 
Poo \ 2il o J 

Equation (J6]) shows that the intermediate state decay 
can be neglected during the transfer time T tr if (Ti + 

r 2 )Tt r (sin20f2/2n o ) < 1- From this expression one 
can see that the intermediate state decay rate is reduced 
by a factor (£1 / S^o) 2 *C 1 in this regime. It also follows 
from Eq. (|6|) that decay from \gi) is not suppressed, so 
that the transfer process has to be faster than this decay. 

Magneto- and photo-association techniques produce 
molecules mostly from ultracold Bose, two-spin compo- 
nent Fermi and mixture of alkali metal atomic gases. 
In traps with high initial atomic density, weakly bound 
Feshbach molecules rapidly decay due to inelastic atom- 
molecule collisions, which were found to be the major lim- 
iting factor of molecule lifetime. Depending on the quan- 
tum statistics of the constituent atoms, the alkali dimers 
show different behavior with respect to inelastic atom- 
molecule and molecule-molecule collisions. Fermionic al- 
kali dimers in the Feshbach state are very stable with 
respect to collisions, especially close to the resonance, 
where the scattering length is large. Lifetimes of the Fes- 
hbach molecules of the order of 1 s have been observed 
experimentally [li|[2(|. In contrast, bosonic and mixed 
dimers experience fast vibrational quenching due to in- 
elastic atom-molecule collisions, even in their Feshbach 
state. An atomic density in a trap is typically in the 
range n at ~ 10 11 — 10 14 cm" 3 , then the Feshbach state 
relaxation rate is in the range I\ ~ 10 1 — 10 4 s" 1 (cal- 
culated from the corresponding inelastic atom-molecule 
collision coefficient k ine i ~ 10" 10 cm 3 s _1 [H HI HI)- 
At the same atomic densities the vibrational relaxation 
rate r 2 of intermediate vibrational states for bosonic 
molecules is in the range r 2 ~ 10 2 — 10 5 s" 1 (calculated 
from kinei ~ 6 • 10~ 10 cm 3 s _1 for 7 Li2 [24| and the same 
range of atomic densities). Inelastic molecule-molecule 
collisional relaxation rates are about two orders of mag- 
nitude smaller due to typically smaller molecular density. 

We next illustrate the technique for a sample seven- 
state bosonic 87 Rb 2 molecular system (see the inset to 
FigHh). In the first step, the Feshbach state can be cou- 
pled to the electronically excited pure long range molecu- 
lar state |0~, v, J = 0), located close to the 5S1/2 + 5P3/2 
dissociation asymptote. For example, following [ll|, the 
v = 31 vibrational level can be chosen 6.87 cm" 1 below 



TABLE I: Possible chainwise transfer path from the Feshbach 
to the ground rovibrational state in the Rb2 molecule. Shown 
also the corresponding transition dipole moments and wave- 
lengths. 

1 v — v' transition ~Dv~vH A 

Debye nm 

1 \Feshbach) ~ |0" , v = 31, J = 0) OA 780.7 

2 |0J,« = 31,J = 0)-X 1 E+(w = 116,J = 0) 0.8 780.4 

3 X 1 T,g(v = 116, J = 0) - J 4 1 E+(w' = 152, J = 1) 0.55 846 

4 A 1 E+(«' = 152,J=1)-X 1 E+(« = 50,J = 0) 0.64 907.4 

5 X 1 Z+(v = 50, J = 0) - A 1 E+(u' = 21, J = 1) 0.53 990 

6 A = 21, J = 1) -X'T.jjv = 0, J = 0) 2.37 856.4 



the dissociation limit. The second STIRAP step can be 
to v = 116 in the ground electronic state. The authors 
of Ref.[ll[ mention that the Franck-Condon factors from 
the excited |0~, v — 31, J = 0) state to the ground state 
vibrational levels down to the X 1 S+(t; = 116) are sim- 
ilar to the second-to-last vibrational state used in the 
STIRAP experiment in [ll| . The ground v = state can 
then be reached in four steps, using e.g. the path given 
in Table I. We note that in Rb 2 the v — state cannot 
be reached from the v = 116 in two steps due to unfa- 
vorable Franck-Condon factors, a minimum of four steps 
is therefore required, resulting in a seven-state system. 
A four-step path from the Feshbach to the v = state 
can be realized in Cs 2 [Hj], and therefore in other alkali 
dimers as well. 

The transitions |ei) - \g 2 ); \g 2 ) - |e 2 ), |e 2 ) - \g 3 ) 
and \g 3 ) — |e 3 ) are coupled by CW laser fields, the 
first transition \gi) — |ei) and the last transition |es) — 
I $4) in the chain are coupled by the fields Cl\ = 
Qrnaxf^ + t anh (t - r/2)/T)/2 and Q 6 = fl^ ax (l - 
tanh (t + r/2)/T)/2, respectively. In the above scheme 
we picked the |es) — | (74) transition with a large transition 
dipole moment, and intermediate transitions coupled by 
CW fields with close and reasonably large moments. In 
this case the Stokes pulse intensity can be minimized, and 
CW fields, provided, e.g., by laser diodes, can have the 
same intensity to optimize the transfer efficiency. The 
wavelengths of the transitions in Table I are covered by 
Ti:Sapphire and diode lasers. To provide the phase co- 
herence between the laser fields required to carry out 
STIRAP, lasers can be phase locked to spectral compo- 
nents of a frequency comb [26| . 

The results of the numerical simulation are given in 
Fig.O We assumed that the CW fields have the same 
amplitude of the electric field £q, resulting in a Rabi fre- 
quency Qq i = D v yi£o/2h for ith transition. The Rabi 
frequencies of STIRAP fields were chosen to satisfy a 
condition that is less than the binding energy of the 
Feshbach molecular state to minimize Raman dissocia- 
tion of weakly bound molecules. The pulse duration T 
and delay r were varied to obtain the maximal transfer 
efficiency. To estimate the decay rate of intermediate vi- 
brational states, the highest atomic density n at ~ 10 14 
cm" 3 available experimentally was used along with the 
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FIG. 2: Results of numerical solution of the density matrix 
equation for a seven-state 87 Rb2 molecular system: a) pump 
Hi (solid line) and Stokes Qa (dashed line) Rabi frequencies; 
b) populations of the intermediate vibrational states; c) pop- 
ulations of the Feshbach p gigi (dashed line) and the ground 
v = state p g3 g 3 (solid line). Parameters used: Ti = 10 4 s _1 , 
r 2 = T 3 = 6 ■ 10 4 s"\ 71 = 8 ■ 10 7 s - \ 72 = 73 = 3 ■ 10 7 

s -l. Q utax = Q rnax = 3 . 1Q 7 g"^ Q. = 1 2 . D v>v l j • 10 8 S" 1 

i=2,3,4,5, T = 1 fis, t = —2 fis for both schemes. Inset in (a) 
shows a model Rb2 system with decays. 



inelastic collision coefficient for intermediate vibrational 
states k ine i ~6-lCr 10 cmV 1 , giving r 2 , 3 = 6-10 4 s~ 1 . A 
Decay rate of the Feshbach state Ti = 10 4 s _1 was used. 
Numerical analysis shows that > 90% of the population 
can be transferred to v — at high initial atomic density 



even in the presence of collisional decay from the intial 
Feshbach state. As can be seen from FigJ^p, the popu- 
lation of the intermediate ground vibrational states does 
not exceed 7% during the transfer process and only for a 
short time, reducing the molecular loss due to collisions 
in this states. 

We can now estimate intensities of CW and pulsed 
fields corresponding to Rabi frequencies used in our cal- 
culations. Taking the peak Rabi frequency of the pump 



and Stokes fields Vt^ x = 



j-pcah 

1,6 



3T0 s , the corresponding in 
3 /8tt = c(n^g*h/D ViV , r, 6 ) 2 
3 W/cm 2 and J 6 ~ 0.1 W/cm 2 



87T, 

for 
the 



tensities are 
resulting in I\ 

CW fields with a Rabi frequency fioi ~ 6 ■ 10 7 
corresponding intensity is ^2,3,4,5 ~ 5 W/cm 2 . 

In summary, we propose a method of vibrational cool- 
ing of ultracold molecules, based on the multistate chain- 
wise STIRAP technique. Molecules which are formed in 
high-lying vibrational states are transfered into a ground 
rovibrational state v — 0, J — using Raman transi- 
tions via several intermediate vibrational states in the 
ground electronic state. Our technique provides 100% vi- 
brational as well as rotational selectivity using selection 
rules A J = 0,±1 for rotational transitions. Numerical 
analysis of the transfer process for a typical bosonic Rb 2 
molecular system in a trap with a high atomic density 
n a t ~ 10 14 cm~ 3 shows that transfer efficiencies ~ 90% 
are possible even in the presence of fast collisional relax- 
ation of the Feshbach molecular state. 

The multistate chainwise STIRAP technique allows 
one to use various transitions, coupled by, e.g., rf fields 
and DC interactions. It can therefore be combined with 
the recently demonstrated resonant association method 
[27j |. Another possibility is to use the magnetic field de- 
pendent DC interchannel coupling between an entrance 
and a closed channel state as a first transition in the 
STIRAP chain [28] followed by optical transitions to the 
ground vibrational state. The chainwise STIRAP can 
be applied to resonant photoassociation as well, then the 
first transition in the STIRAP chain will couple the con- 
tinuum states to a high energy vibrational state in the 
ground electronic state [29(. 
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